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Introduction 

A precise geopotentisl aodel has becooe increasingly inportant as more 
satellites, particularly those designed for geodetic and altlnetry purposes, 
require accurate orbital tracking to produce meaningful resulta from their 
orbital data. The geoid model, the static part of the geopotential field. has 
been rapidly improved using various data from satellite tracking and surface 
gravity iu terms of the order and degree of spherical harmonic expansion. Most 
recent models may include NASA's Goddar.1 Space Plight Center geoid based on the 
work of Marsh and Vincent (1974), GEM 9 and 10 by Lerch et al (1977), and one re- 
ported by Gaposchkin and Mendes (1977) , the last two models having the degree 
anu order lO. 

Most importantly, the geoid model cannot be significantly improved until we 
know the distribution of the open ocean tides. Since tides are regarded as noise 
in deriving a geoid model, tht RMS error amplitude of the resultant geoid can 
not be theoretically smaller than the RMS amplitude of tides. Those geoid 
models having residual amplitudes less than a few meters may have already reached 
this barrier. To break this barrier, we are inevitably led to consider a tech- 
nique which is capable of simultaneously deriving both geoid and tides 
from satellite altimeter data. 

Satellite altimeter data consist of ocean elevation measurements with 
respect to a reference spheroid along subsatellite tracks distributed more or 
less randomly both in space and time. The mapping of the tides using these 
data is considerably more difficult both mathematically and computationally 
than that required for a static geoid and requires a far greater amount of 
altimeter data with higher accuracy. This is not only because of the dynamic 
nature of the tides, but also their relatively sDiall amplitudes, generally less 
than one meter worldwide. 

Zetler and Maul (1971) performed a study using simulated satellite alti- 
meter data on the tidal analysis at a fixed location. Tl\ey assumed that the 
tidal height at a given time would be uniform ever a 5* square area and that they 


- 1 - 


coulc! use the altimeter data whenei^er the satellite passed over the area. A 
randomly sampled time series thus generated may he solved for tides hv the 
least-squares harvonlc method. Thev shoved that It was possible to retrieve 
ouijor tidal constituents from the simulated data In the presence of random 
noise which was larger than the signal. 

Maul and Yanawsy (1977) applied this method to nKOS-3 data obtained within 
a 5* X 5* squa'.e centered at ?0*N and 70*W (near Bermuda Island). Thev 
regarded the data falling In this square as a time-series of a single location 
and applied a least-square harmonic analysis method which Is regularlv used 
for conventional tidal analysis. The results they obtained showed that the 
tidal amplitudes from the GF.OS-3 data were an order of magnitude higher than 
those obtained from the MODE deep-sea tide gauge, and their solution did not 
converge when the size of the square was changed. Therefore, thev concluded 
that It was not possible to retrieve tides from GEOS-l data. 

Won et al (1977) developed a similar least-square harmonic analysis method 
which Is designed to simultaneously extract both the oceanic geold and tides 
for a surface (rather than for a single point) using satellite altimeter data. 
In their case, both the geold and tides within a region are represented bv a 
set of two-dimensional functions. Ry applying this method to sets of simulated 
data In the northeastern Pacific Ocean, they showed that tides whose amplitude 
may be as small as 10 cm nwy be retrieved from the simulated data even when the 
data .are contaminated with meter random noise. 

Both studies described above assume that the error In the altimeter data 
Is random and has no orbital drift. However, the most severe dltficultv 
In using altimeter data for the tidal analysis stems from orbital bias error. 

Marsh et al (1976) Investigated the problem with the Skylab altimeter data, 
and found the error is position-dependent and Is of long wavelength nature 

(*100®) with amplitudes on the order of several mi>ters. Tl>e orbital bias 
error of GEOS-3 has been steadily Improved since Its launch and stahili;'etl 
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Co !••• CliAn • Mtor In the Calibration Area (U. R. Stanlay, paraonal 


conunlcaClon, 1977). 


DaCeralnatlon of ocean Cldea for Cha anClra northam Atlantic may require 


aevaral hundred patha for a duration of a year or more, baaed on the alnulated 


atudy by Won at al (1977). Since available data on hand are limited, we 


preaent In chla paper the reaults of the oceanic geold and four major tides. 


^2* ^1* ^2* ^l* along tvo linear patha of CEOS-3 as shown In 


Figure 1. We shall briefly dlscuas the analytic method, results of a simulated 


study Co test the method, and finally the results of actual CEOS-3 data. 


Analysis of Linear Satellite Data for Oceanic Ceold and Tides 


The conventional tidal analysis technique requires a long time 'er^es 


obtained at a fixed point to obtain Its tidal spectrum. Since Che al Imeter 


oats are obtained more or less randomly In time snd space, this method Is not 


applicable. The following method follows closely Chat by Won et al (1977) 


originally developed for a two dimensional tidal analysis. 


Let us consider H(x,C) the ocean height with respect to a reference 


spheroid at a distance x away from an arbitrary origin along the ground path 


of a satellite. Neglecting transient sea surface changes, we write 


H(x,t) - C (x) + L A. (x) c« {w.t + K (x)} 


where C^j(x): constant ocean height, l.e., ocean geold. 


A^(x), Kj(x): amplitude and phase of 1-th tidal component, 

u)j^: frequency of 1-th tidal Imponent , and 


t; time after the start of measurement 


Since Kj^(x) Is continuous only In the Interval of 0 to 2tt, we use Instead 


C. (x) ) (cos X (x) 

f- 

S^(x)\ ^In Xj(x) 
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Figure 1. Locations of the two strips of data In the GEOS-3 calibration 
area used for the analysis. The origin for counting linear 
distance la shown as a solid circle for each strip. The width 
of the strips is approximately 265 Km at 30*N latitude. 




anti l«t and S^, r««p«tct Ivr ly, ha rapraarntad bv a polvnoalal aarlaa auch that 
whara T (x) la anv alaaantary function of ordar a In x. 


For a glvan aat of N altlnatar data I’fx^.t^); j • 1,2 
wlah to find tha haat fitting H(x,t) ao that tha quantity 


. . N, wa 


S - r lF(x ,t ) - H(x ,t )1^ (4) 

j.l ' J ' ' 

ha alnlaltad with raapoct to ami a^. 

Whan tha gaold and tldaa, raapact Ivaly, ara approxlmatad hv and ordar 
funrtlona and L Indapandant tidal rimponanta ara aought, tha aquation (41 
ganarataa a poalt Iva-daf Inl ta avamatrlc matrix of ordar + 2N^L. 

Thara ara two aarloua aourcaa of maaauromont arror (Brvjw Ct aJ l'*7h> 
which can complataly Invalldata tha achama praaantad hara. Tha flrat aourca 
la, of couraa, ihr ramliWB moHaurrxiant arror of Inatrumaut. Howavar, aa ahowo 
by Zatlar and Haul (1171) , at_ aj|_ (117 7) , Maul (1177) , and In thla rap»>rt , 
thla arror ran ha ovarcoma aa long aa It la random. 

Tha aacond arror aourra, much XK'ra aarloua and aaamlnglv unraaolvahla, la 


that cauaad hv tha unrartalnty of aatalllta orhlt, fraquantlv rafarrad to aa 
orbital hlaa arror. Kot ax.impla, tlia lasulta radurad tiom tha iIKOS- 1 .iltlmatiM 


data In the Calibration Area ahi^ that tha aaa aurfara halghta along a glvan 
ground path, whan roraparad with a nx'dal gaotd, deviate aa much aa 20 matara. 

In addition, linear drift la t>ptcallv about ) itu»tara par 1000 km along a 
ground path. Tlta errors ara nvilnlv attributed to uncart alnf las In gaopotantlal 
uK'dala used for .atalllta tracking. 

Tha critical problem In tidal analvala la that these orbital bias errors 
will be aliased Into anv tidal spactrimi. Tha problem can be more serious If 
tha orbital parlixl of a satellite Is In resonance with anv tidal period. I'ne 
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poHslhle aethod la Incorporating the orbital blaa error Into equation (4) by aa- 
aualng that the hlaa la repreaented by a polvnoalal for each path. In other 
worda. Instead of (4), we alnlalze 

“ H ^ « "a*?’" 

where the subscript p denotes quantities In the p-th path only. 
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R esults of Analysis of Slaulated Data 

To construct a realistic ocean surface as a function of space and time, we 
used a satellite geold to the degree and order 20 bv Rapp (1474) and aaplltude 
and phase aaps of M 2 and 0^ tides by Hrp.O fit a 1 (1467) In the northeastern 
Pacific Ocean. We arbitrarily chose one path similar to that of GF.OS-1 and 
reconstructed the sea surface height data along this path. The geold undula- 
tion along this path amounts to about 30 m, while the tides ate less than Im 
everywhere. For the given GEOS-3 orbital period (102 minutes) and Inclination 
(65*), the data are collected whenever the mathematical satellite passes over 
within 2* In longitude on uitlier side ot t lie equator crossing of the 
designated path. Sampling time was chosen to e 3.2 seconds, (approxlmatelv 22 km 
in distance) similar to that of GKOS-3. One data set thus generated Is comprised 
of 2000 altimeter data points from 16 paths, each path about 3000 km long, 
covering 24.2 days. 

A random series having a two meter peak-to-peak arni^lltude range is then 
generated in order to artificially contaminate the simulated data. Tlie white 
noise series has a mi*an of near zero and an RMS amplitude of 58.4 cm. Tlie 
series is then piecewise added to each path on a polat-by-point basis. 

When the noise was not added to the data, the reproduced geold and tidal 
profiles were almost Identical to the prescribed ones as shown In solid lines 
In Figure 2. The figure also sliows results wlien the data were dubbed with 1 
meter (dashed lines) and 2 meter (broken lines) noise. 
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In all cases, the reproduction of the geold Is assured. While tidal 
signals are nnstly well below the RMS noise amplitude, their recovery Is 
reasonably good. Several factors contribute to the successful recovery: 

(1) the added noise Is Hausslsn, (2) no artificial drift within each path 
(such as the orbital bias error) Is allowed, and (3) tidal frequencies are 
precisely known. 


Description of the CKOS-3 Data Used for the Analys is 

Two sets of GEOS-3 altimeter data confined In two narrow strips (Figure 1) 
are chosen for ',)resent analysis. One strip, starting at Newfoundland and ending 
at Cuba, consists of 28 south-going paths whose equator crossings are between 
99.64*E and 102.46*F.. The width of tne strip Is about 284 km at the equator 
and narrows somewhat toward higher latitude. The other strip, starting from 
near Puerto Rico and ending at the North Carolina coast, consists of 26 north- 
going paths whose equator crossings are between 53.54*W and 56.36*W. The data 
duration for both paths Is about 17 months, start Ing April, 1*^75, and tepresent 
about 701 of tile total available data within these two strips during that 
period. 

Tlie first two graphs of Figures 3 and 4 show the original altimeter sea- 
surface height data along with geold heights computed from the Goddard Space 
Flight Center Ma'sh-Vlnrent geopotential model. It Is evident that the 
orbital bias error Is large, compared with the computed geoldal height 
frequently exceeding 10 ra or more. 

The Initial data preparation for each strip proceeds as follows: 

Step 1. Select a hypothetical reference path through the middle of the strip 
and project each data point onto the reference path along a line on 
which the dlntance between the data point and the reference path Is 
mini mum. 
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Step 2. Assign sn arbitrary origin ou the reference path a. compute the great 
circle distance for each data point %rith respect to thia origin. This 
step assluilates the original altimeter data into a single aet of 
linear data aa a function of distance. 

Step 3. Remove bad data poirts. 

Step 4. For data which fall in a specified interval distance, compute and 

remove a beat-fitting linear function for each path. The results are 
shown in the third graphs of Figures 3 and 4. 

It should be pointed out tiiat during the last step any geoldal or tidal 
fluctuation whose wavelengths are longer chan the spatial length is irretriev- 
ably lost. Furtherinure, any nonlinear orbit drift is not removed and will 
affect the ''ubsequent analysis. Wliile any geoidal fluctuation %diose wavelength 
is greater than about 1200 k.m is supposedly included in the Marsh-VlncenC geold, 
simi''ar argument may not hold for the tidal fluctuation. 

The characteristics of the GEOS-3 data used for this analysis are 
summarized In Table 1. The locations of the reference paths are also shown 
in Table 1 as well as in Figure 1. The final data set resulting from this 
preliminary preparation is then used for the least-square space-time har- 
monic analysis as described in the previous section in order to simultaneously 


obtain the tlder and the residual ocean geold. 


TABLF. 1. CharacttfrlHtlCH of c;K 0S-3 Data Uaed for 
Analyala and tho Local Iona of Reference Patha 


SOUTH-GO INC STRIP NORTH-1'.OING STRIP 


a. Number of patha uaed 

h. Duration 

c. Total uaable data 
point a (coverage) 

d. Riinge of equator 
croaalng anglea 

e. Strip width 

equator 
30* latitude 

f. M**an offael from 
geo Id* 

g. ^U•an a lope with 
reapect to geo Id 

h. Kquator croaalng 
of reference path 

1. Origin for linear 
dlatance along the 
reference path 


Id 

Miiy Id75-Aug 197b 
4591 

(New found land-Guba) 
99.b4*K-102.46*K 


284km 

265km 

9.H5m(SD'*^-4.23m) 


-4. 5 1 m/ 1000km 
(SD+-1 .77m/lOOOkm) 

100.«700'K. 


55.2809*W 

49.2387’N 

(Newfoundland) 


26 

Apr 1975-Aug 1976 
2929 

(Puerto Rlco-N .Carol Ina) 
53.54*W-56.36*W 


284km 

265km 

-5.54m(Sl/-3.74m) 


1 . 00m/ 1 000km 
(SlV‘’-l.19kra/ 1000km) 

55.0100’W 


62.1475*W 
12.8496“N 
(N. of Trinidad) 


* Miirsh-Vlncent lUiddard Space Flight Center Ceold 

+ Stamlard deviation 
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lUiiult of GEOS-3 Data Analysia 

Th« present anelysls teclinlque hen e few erbitrery yet Inportent 
variebles which may produce different solutions for e given set of date. 
These variables include (1) nuaber of tidal constituents sought (In addition 
to a geold), (2) order of polynomial function representing i geold as 
well as tides, and (3) spatial length of data. While the first two 
variables are attributed to the usual convergence problen resulting froa 
Incompleteness of both representation function and data, the third variable 
Is due to the removal of the orbital drift error as described In Step 4 
In the previous section. 

Since these variables aay significantly affect the convergence of the 
solution, we allowed then to vary within limited. rangmUk fi(>eclf Ically, 
analyses were performed for the following range of each variable: 

1. Number of tidal constituents - 2 groups 

1st Group: Geold, and 0^ tides 

2nd Group: Geold, 0^, S^, and tides 

2. Order of polynomial - 4 groups 

5th, 6th, 7th, and 8th order Chevychmy polynomials for both 

geold and tides 

3. Spatial length of data 

a) 4 groups for the south-going paths: 300 iCm . 1,600 Km; 

300 Km ' 2,400 Km; 2,200 Km - 3,600 Km; and 300 Km > 3500 Km 

b) 1 group for the north-going paths; 1,000 Km . 2,900 Km 

In total, 32 different analyses were performed for the south-going paths, 
and 8 for the north-going paths. An origin time for the entire data was 
arbitrarily chosen to Le 00:00:00 GMT on January 1, 1975 so that all tidal 
phases are referenced to this time. 

-13- 
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Incorporation of the orbital drift error into the norm.il equation (5) 
was tested for both north- and sout'i-ftoinK data sets. When the orbital 
drift for each path is expressed by a linear equation in distance, the 
process K^nerates additional 2N^ unknowns where is the total number of 
paths used for analysis. The data preparation tor this test, of course, 
did not include Step 4 described in the previous section. Test results 
showed that the coefficients of the linear equation for orbital drift thus 
obtained were almost identical to those obtained from Step 4. However, the 
increased number of norioal equations caused numerical instability for high- 
order (7th and 8th order) polynomial representations of the geold and tides 
due to the limited amount of input data. Therefore, this process was not 
pursued any further. 

It is Impractical to present all 40 different sets of solutions in 
this report; we present only a set of four different solutions for each 
strip in figures 5 and ft. Other solutions show similar tidal characteristics 
without significant di ♦’ferences. A statistical summary for all 40 solutions 
is shown in Table 2. 

For a given spatial length of data it was found that the geold solutions 
agree within 30 cm or less for all orders of polynomial and more or less 
Independent vif the number of tidal constituents. Change in the spatial data 
lengt 1 caused only "tilting" of the geold, mainly because the length aifects 
the linear fitting (Step 4) during the data preparation in removing the 
orbital drift errors. 

The residual geolds obtained from the analysis lor both south- and north- 
going paths are plotted on the fourth graph of Figures 3 and 4. Eacli solution 
is compared with the simple arithmetic mean of the residial altimeter data. 

In both c.ises, the agreement appears to be re.isonable. 
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Figure S. South-going amplitude and phase profiles for M2, Oj , S2 *nd I, tides obtained fro« altlvcter 
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Figure 6. North-going amplitude and phase profiles for M 2 , 0^^, S 2 

tides. Data from MODE gauge(Flgure 1) are shovm In solid circles, 
additional explanation refer to Figure 5. 
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TABLE 2. RMS Error Aaplltudefl(ca) of che Residual Alclaetcr 
Data after Reaovlog Geold and Tides for AO Different Analyses 
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Wlu*n the reHldual xeoid than obtained 1 h added to the original M.irHh- 
Vlncent geoid, we obtain an "Improved" geoid eontalning short w.ive length 
features which are not Included In the original geoid model. The "Improved" 
geolds are shown on the fifth graph of figures 3 and 4, and compared with the 
Marsh-V Incent ones. It Is interesting to note somewli.it better correlation of 
the "Improved" geoid with bathymetric profiles. For the north-going paths, 
the geoldal correction amounts to as much as b m near the Puerto Rican Trench, 
however, this large deviation nuiy be partially due to a relatively short 
spatial data lengtii along these paths. 

We note from Figures 5 .ltd b th.it the four solutions for tidal amplitudes 
agree within about 10 cm in most cases except near nuirgins due to familiar 
edge effect. Agreement on phases is somewh.it poor although consistent trends 
are noticeable In .ill cases. 

Figure b .ilso slu>ws the tidal amplitudes obtalneil from the MOPK deep-sea 
gauge (feller ^ aj , ld7!>). Wlille .ill t id.il amplitudes .ire within .in order 
of magnitude not exceeding 40 cm, the agreement with the MODF. d.ita is r.ither 
poor particularly for M,, tide. It Is disturbing ti> note th.it S., amplitude is 
bigger th.in that of M., in figures S and b. l^onsidering t lu> tid.il inducing 
force, we expect M, to h.ive gre.iti-r .implitude th.in unless there exists an 
unusual spectral feedh.ick between the two i-onst i t uent s in the dat.i area. It 
is miire likely ilue to the Insuftlcient ipi.intlty .ind .n'cur.icv ot .iltlmeter d.it.i 

l)n till* other li.ind, .inv griMiiul truth d.it.i from di'ep se.i iiuMsurement .mu 
be Incorporat Oil in the norm.il equations using the l.agr.aigl.in multiplier tech- 
nique iWo^ £t^ 1477). This technique w.is not applied in this .in.ilysis 

m.ilnly because the presiuit |nirpose i.s to compare the two dilleriuit solutions. 


It is anticipated that aa more data are accumulated from GEOS-3 and 
future SEASAT aerlea aatellltea with Improved tracking technique, the dis- 
crepancy will resolve eventually. The main emphasis here must be that the 
analytic method presented here does produce a consistently converging solution, 
within expectation. In spite of limited quantity and quality of the altimeter 
data. 

Conclusions 

1. Happing the ocean tides and geoid using satellite altimeter data 
is indeed possible In a large regional scale provided that a sufficient amount 
of data and reasonable altimeter accuracy are attainable. 

From the viewpoint of economy and available data density, the satellite 
altimeter provides the most abundant and direct measurements. The satellite 
altimeter data can best be utilized in obtaining a geopotential model when 
the data are reduced simultaneously for a geoid and for a set of tidal surfaces. 
Large residual errors which are present In available gravity models are mainly 
attributed to: 1) exclusion of the tidal perturbation and, 2) ineffectiveness 

of the spherical harmonic expansion technique in describing localized gravity 
anomalies. These two factors, in addition to the measurement errors, essen- 
tially dictate the upper limit of the residuals of a gravity model, which may 
be on the order of a few meters or less. Some of the recent gravity models are 
already approaching this limit. 

Considering the large number of tidal constituents, the regional treatment 
of the altimeter data appears to be the best approach in deriving botli the 
static and the dynamic geopotential models. The entire ocean area nuty be 
divided into about a dozen partially overlapping regions. For each region, 
geoid and tide models will be derived from the altimeter data by specifying 
a functional surface for each model whose resolution can be much higher than 
those att.iinable by the spherical harmonic technique. 
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2 . Th« ocennlc can Sr laprovril alunlf leant ly throuith ( l^lal 

analyala. In a atrlct arnar, the ocean gei>ld obtained from the altlneter 
data la a t In*- Invar lant ocean aurface topography with respect to the 
spheroid. Therefore the ocean geold thus derived, even after tidal 
corrections, laav still be contaalnated bv localised elevations created 
bv steadv-state geostrophlc current svsteaa. If the current svsten 
meanders fast enough as compared with the duration of the altimeter data, 
we may hope that Its long-term effect on the geold would not be sig- 
nificant. In addition, since the geostrophlc currents occupy only a 
small portion of the ocean, the overall residual amplitude of a geold can 
be .educi'd l»’ slgnl I leant 1 v less than one meter when tides are fullv 
accoimted for. 

d. three the altimeter data are corrected for the geold and tides, the 
resldikil altlns'ter data can be used to Investigate tr.rnslent sea surf.ice 
variations wlilcli nuy be associated with winds and storins. Tlte residual 
data uwiv also be used to Identify exlretm'lv locallr.ed geoldal 
features caused bv snwill sea tm'iuits and narrow trendies. 
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